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SEMIEMPIRICAL  POTENTIAL  ENERGY  SURFACES  OF  (ArH20)+  SYSTEM 
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Abstract 

In  order  to  study  the  exothermic  charge  transfer  reaction  Ar“^  +  H2O  Ar  + 
(H20)'^  the  semiempirical  calculation  of  the  (Arll2  0)“^  potential  energy  surfaces 
(PESs)  is  performed.  These  PESs  are  used  to  calculate  the  vibrational  modes  of  the 
molecular  species  involved,  in  particular  the  highly  excited  vibronic  bending  modes 
of  the  reaction  product  (H20)'^.  It  was  found  that  the  bending  modes  are  affected 
by  the  stretching-bending  coupling.  It  is  shown  that  at  suprathermal  collision  ener¬ 
gies  the  electrostatic  approximation  is  not  valid  anymore  as  the  electron  hop  occurs 
mainly  in  the  strong  interaction  region.  According  to  the  results  of  the  full  scale  cal¬ 
culation  the  reagents  Ar“^+H2  0  vibrational  ground  state  is  in  near  resonance  with  the 
excited  bending  (0,9,0)  state  of  the  products  Ar+{H20)‘^  along  the  whole  interaction 
region.  In  a  derivation  based  on  simplifying  cissumptions  (R.A.Dressler,  J.A. Gardner, 
R.H.Salter,  and  E.Murad,  J.Chem.Phys.  1062  (1992))  the  PESs  were  found  to  be 
significantly  different  in  particularly  the  Ar'^+H20  vibrational  ground  state  crosses  the 
Ar+(H20)“^  bending  manifold  in  the  strong  interaction  region. 
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1.  Introduction 


The  purpose  of  this  work  is  to  calculate  the  (ArH20)"*"  potential  energy  surfaces 
(PES)  which  describe  the  charge  transfer  process  for  the 

Ar+  +  H2O  ^  Ar  +  (H2  0)+  (l) 

reaction.^  The  charge  transfer  process  between  rare  gas  atoms  and  molecules  have  re¬ 
ceived  during  the  years  considerable  attention.  However  the  theoretical  studies  in  this 
field  have  mainly  emphasized  diatomic  molecules  and  consequently  systems  such  as 
(ArH2)+,2-8  (ArN2)+,^-i°,  (ArCO)+,9  (HeH2)+ ®  (NeH2)+,®  and  (HeN2)+  were 
considered.  No  such  studies  are  known  to  us  for  systems  with  triatomic  molecules  ex¬ 
cept  for  (HeHs)^  and  (ArH3)+  Our  study  is,  probably,  the  first  one  dealing  with 
the  interaction  of  a  rare  gas  ion  with  a  triatomic  molecule,  such  as  H2O. 

The  interaction  between  a  rare  gas  ion  (atom)  and  a  polar  molecule  (molecular 
ion)  is  easily  calculated  in  the  asymptotic  region  where  the  interactions  are  determined 
by  electrostatic  forces  such  as  charge-dipole  forces  in  the  case  of  the  Ar+-t-M  system 
and  polarization  forces  in  the  case  of  the  Ar-f-M"^  system.  However  the  polarization 
potential  decreases  very  fast  as  R  increases  (namely  like  l/R^)  and  so  the  charge-dipole 
interaction  which  decreases  with  R  relatively  slowly  (as  l/i?^)  becomes  the  dominant 
interaction  in  this  region.  This  fact  was  used  by  Dressier  et  al.^  to  construct  their 
approximation  for  the  (ArH20)+  which  assumes  the  charge-dipole  interaction  to  be 
valid  for  the  whole  range.  However  this  approximation  is  not  valid  at  relatively  small 
distances  where  interactions  resulting  from  orbital  overlap  and  charge  delocalization 
become  strong.  The  calculations  of  these  interactions  have  to  be  performed  either  by 
an  ah  imtio  treatment  or  by  semiempirical  methods.  In  this  work  such  a  semiempirical 
treatment  is  presented. 

The  PES  calculations  of  reactive  systems  are  carried  out  mostly  by  employing 
diatomics-in-molecules  (DIM)  method.^-^-i®  The  main  advantage  of  this  method  is  that 
it  offers  the  possibility  of  expressing  the  energy  of  a  polyatomic  system  in  terms  of 
diatomic  potentials.  The  application  of  the  DIM  method  to  ionic  systems  such  as 
(ArH20|'^  presents,  however,  some  diflSculties,  mainly  because  of  the  charge  delocaliza¬ 
tion  and  electrostatic  interactions.  These  difficulties  are  successfully  removed  employing 
the  diatomics-in-ionic-systems  (DIIS)  method  which,  like  the  DIM  method,  applies  di¬ 
atomic  potentials  to  produce  polyatomic  PES.^^'^o  y^^ereas  within  the  DIM 
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one  employs  adiabatic  potentials  only,  within  the  DIIS  one  incorporates  both  adiabatic 
as  well  as  diabatic  diatomic  potentials.  This  modification  enabl^the  DIIS  method 
to  treat  simple  systems  with  charge  delocalization,  such  as  ionic  rare  gas  clusters^  or 
excited  rare  gas  —  halogen  molecules  Xe^Cl’".^^ 

The  charge  transfer  reaction  is  described  as  an  nonadiabatic  electronic  transition 
and  consequently  at  least  two  PESs  are  involved  in  the  process.  The  charge  transfer 
reaction  (1)  leads  to  the  formation  of  the  (1120)“^  molecular  ion  in  the  excited  electronic 
state  A^Ai.  This  is  known  to  be  a  highly  exothermic  reaction  and  consequently  the 
(H20)+  ion  is  expected  to  be  formed  in  highly  excited  vibronic  bending  states.^  It 
follows  that  in  order  to  understand  the  charge  transfer  process  in  the  (ArH20)‘*‘system 
we  need  to  calculate  not  only  the  (ArH20)'*'  PESs  but  also  the  vibrational  modes  of 
the  neutral  and  ionic  species  along  the  reaction  coordinate. 


2.  Semiempirical  treatment  of  the  (ArH20)+  PES 

As  was  mentioned  earlier  the  wave  functions  within  the  DIIS  are  expressed  in 
terms  of  adiabatic  and  diabatic  diatomic  wave  functions.  The  diabatic  wave  functions 
describe  electronic  configurations  with  localized  charges,  such  as  Ar+H  and  ArH+  but 
not  (ArH)+  which  are  presented  in  terms  of  the  adiabatic  functions. 

In  accordance  with  the  DIIS  method  the  polyatomic  basis  functions  $  are  presented 
in  terms  of  antisymmetrized  products  of  atomic  orbitals  (AO)  Xi  where  *  is  the  atom 
number.  The  (ArH20)+  atoms  are  numbered  as  follows:  the  0-atom  is  i  =  1,  the 
and  Hf)  atoms  are  i  =  1  and  3,  and  the  Ar  atom  is  i  —  4.  The  s-AO  is  denoted  by 
the  atom  number  only  whereas  the  p-AO  is  denoted  by  both  the  atom  number  and  the 
orientation  index.  For  example,  the  Is  AO  of  the  atom  is  Xz,  but  the  2p  AOs  of 
the  O  atom  are  xii,  Xiy.  Xiz-  We  also  distinguish  between  the  two  spin  states  x  (for 
s  =  i)  and  X  (for  s  =  -i) 

As  it  is  usually  accepted  in  semiempirical  methods,  the  O  atom  is  described  by  the 
outer  2p-AOs  ignoring  the  inner  Is  and  2s  electrons.  In  analogy  with  O  the  Ar  atom  is 
described  m  terms  of  the  outer  3p-AOs.  In  addition  we  do  an  important  simplification 
in  which  the  P-symmetry  of  the  Ar+  is  ignored  and  active  valence  AO  3p^,  for  example, 
is  replaced  by  an  s-AO.  Consequently  the  Ar”*"  and  Ar  are  presented  as  X4  aJid  X4X4> 
respectively. 
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While  describing  the  polyatomic  basis  functions  we  omit  the  symbol  of  antisym- 
metrization  operator.  The  basis  functions  are  built  in  such  a  way  that  in  the  asymptotic 
case  of  only  two  interacting  atoms  the  interaction  potentials  coincide  with  the  corre¬ 
sponding  diatomic  potentials.  The  (ArH20)'^  total  spin  ]s  S  =  ~  and  its  projection  is 
assumed  to  be  5^:  =  The  coordinate  system  we  are  using  is  as  follows  (Fig.l):  the 
coordinate  center  is  located  at  the  atom  O,  the  z-axis  is  directed  toward  the  Ha  atom, 
and  the  x-axis  lies  in  the  HOH  plane.  The  location  of  the  Ar  atom  is  arbitrary. 

Before  treating  the  (ArH20)”^  system  as  a  whole  we  consider  the  two  molecular 
fragments  of  this  system,  namely  H2O  and  (H2  0)'^.  The  water  molecule  was  treated 
already  by  the  DIM  method^^’^^-^^  but  we  can  not  use  this  representation  as  in  our 
treatment  H2O  (like  (H20)"^)  is  not  an  isolated  molecule  but  a  fragment  of  a  more 
complicated  system.  This  imposes  restrictions  on  the  wave  function  presentation  and 
the  calculation  procedure.  Treating  semiempirically  H2O  and  (H20)+  we  did  the  fol¬ 
lowing:  (a)  a  single  set  of  basis  functions  and  diatomic  potentials  were  used  to  describe 
simultaneously  both  molecular  fragments,  (b)  the  calculated  parameters  of  the  grotmd 
{X  ^Ai)  H2O  and  the  ground  {X  ^Bi)  and  excited  {A  ^Ai)  states  of  (H20)+  fit  as 
much  as  possible  the  corresponding  empirical  values,  including  the  energy  gaps  of  these 
states,  (c)  the  H2O  and  (H20)'*"  basis  functions  are  chosen  in  such  a  way  that  they  can 
also  describe  the  interactions  with  Ar  and  Ar+,  (d)  the  number  of  basis  functions  is 
kept  as  small  as  possible. 


2.1  Ground  state  H2O 

In  the  usual  language  of  molecular  orbital  theory  the  H2O  ground  state  configura¬ 
tion  is: 

(Ia)'^(2a)2(l6^^)2(3c)2(i6j^)2  (2) 

where  a  and  b  are  symmetrical  and  anti-symmetrical  orbitals,  respectively.  We  describe 
the  electronic  structure  of  the  X  ^Ai  state  by  the  following  basis  functions: 


=  XlxXlyXlj,Xl2X2X3. 

(3) 

^2  =  XlxXlyXlj,Xl^Xl^X2» 

(4) 

^3  =  XlxXlxXlyXlyXlzX2^ 

(40 

^4  =  XlxXlyXiyXlzXizXz^ 

(5) 
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=  XlxXixXlyXiyXlzXz^  (5') 

The  first  function  correlates  with  the  asymptotic  0(3p)HaHi,  configuration.  The  func¬ 
tions  i  2  5  describe  the  configurations  with  the  electron  transfer  from  one  of 

the  H  atoms  to  the  O  atom  and  correlate  with  the  asymptotic  0~(2p)HaH^  (^2,^3) 
or  0-(2p)H+H6  ($4,  $5).  In  our  representation  the  and  H+Hfe  interactions  are 

described  as  diabatic  potentials,  in  contrast  to  the  usual  DIM  approach  where  the  adi¬ 
abatic  potentials  are  used.  Also  some  of  the  0-H  potentials  will  be  assumed  to  be 
diabatic  ones. 

The  diagonal  matrix  elements  of  the  basis  functions  (3)-(5)  are  as  follows: 

Hl,l  =  +  UoH^C^)  -h  (6) 

Hi, 2  =  +  sm^eUo-jj+Ci:)  +  cos^  e  Uo-jj^Cn)  +  Uo,  (7) 

Hz, 3  =  Wo-hA^^)  +  +  sin2|917o_^+(2n)  +  Uq,  {7') 

H4,4  =  +  sin^  0  +  cos^^  +  Uq,  (8) 

H5,s  =  +  cos^e  Wo- +  sin^^  +  Uo.  (8') 

Uo  =  —  Ih  +  Ao,  (9) 

where  0  is  the  HOH  angle,  the  W’s  and  the  U’s  axe  the  adiabatic  and  the  diabatic 
(with  fixed  charge  location)  potentials  respectively  (the  diatomic  states  are  shown  in 
parentheses) ,  I^-  is  the  ionization  potential  of  H  and  Ao  is  the  electronic  affinity  of  O. 

The  diabatic  potentials  U  are  derived  from  the  appropriate  adiabatic  potentials  in 
the  way  suggested  within  the  DIIS  method.^^“20  consider  first  the  homonuclear 

fragment  .  In  the  framework  of  the  ordinary  approach  the  adiabatic  states  of  this 
ion  are  Eg  (attractive)  and  E^  (repulsive)  which  are  the  solutions  of  the  equation 

U-E  V 

V  U-E  “  °  (10) 

Here  U  is  the  diabatic  H+H  potential  (the  state  symbol  E  is  omitted)  and  V  is  the 
diabatic  coupling  term.  Solving  the  inverse  problem,  i.e  expressing  U  and  V  in  terms 
of  the  two  eigenvalues  E  =  Wg,Wu,  one  obtains 

V  =  \{W,+W^),  (11) 
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V  =  1  (IV,  -  ir„), 


In  the  case  of  the  heteronuclear  fragment  OH  there  are  two  (different)  diabatic 
potentials,  namely  Uqh  and  Uq-  ii+  ,  which  fulfill  the  equation 


UoH  —  E  Vqh 
VoH  Uo-H+  —  E 


(12) 


The  corresponding  eigenvalues  are  the  adiabatic  potentials  of  the  ground  state 
Woh{X^^)  which  correlates  asymptotically  with  0(®P)H(^5)  and  the  excited  state 
Woh{2^TI)  which  correlates  asymptotically  with  0“(^P)H+.  The  exchange  term  Vqh 
is  assumed  to  be  a  fitting  function  which  will  be  determined  later.  Solving  Eq.(12) 
yields  for  the  diabatic  potentials  Uqh  and  Uo-h+  the  expressions: 

UoH  =  ^[Woh{X^TI)  +  l^0H(2'*n)  -  (13) 

Uo-H*  =  I  [WoniX^n]  +  WoH{2^n)  +  (13') 

where 

=  \l[WoH{2m)  -  WoH{Xm)]^  -  (14) 

The  exchange  term  Vqh  has  to  be  small  enough  to  provide  a  real  value  for  Wd. 


Using  the  basis  functions  (3)-(5)  one  obtains  the  following  non-zero  oflF-diagonal 
matrix  elements: 


P'1,2  =  COS^  eVoH^, 

(15) 

Pi, 3  =  sin^  OVoHb, 

(15') 

El, 4  =  VoHai 

(15") 

■^2,3 

=  sine  cose  [Uo-H+Ci:)  +  Uo_^+(2n)] 

6  0 

(16) 

H4,5 

=  sind  cos0[H'o-Hi(^2)  -h  t^o-Hfc(^n)] 

(16') 

II 

to 

II 

(17) 

where  Vqh  and  Vff+  are  the  exchange  terms  of  Eqs.(12)  and  (10),  respectively. 
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2.2  Ground  state  (H20)+ 


The  (H20)‘'‘  ground  state  configuration  is  (for  notations  see  Eq.(2)): 

(Ia)2(2a)2(l6^,)2(3a)2(l6j,)i  ,  (18) 

We  describe  the  electronic  structure  of  the  X  state  by  following  basis  functions: 


—  XlxXlyX\zX\zX2’i 

(19) 

$2 

—  X\xX\xXlyXlzX2’i 

(19') 

CO 

=  XlxXlyXlzXlzXz^ 

(20) 

$4 

—  XlxXixXiyXizXzt 

(20') 

$5 

=  XlxXlxXlyXlzXlz^ 

(21) 

The  pairs  of  functions  ($i,  $2)  and  ($3,  $4)  correlate  with  the  asymptotic  0(®P)HaH^ 
and  0(®P)H+H6  configurations,  respectively,  and  the  function  $s  correlates  with  the 
asymptotic  0“(^P)H+H^  configuration.  The  diagonal  matrix  elements  for  the  basis 


functions  (19)-(21)  are  as  follows: 

Hi, I  =  WohA^^~)  +  cos2tfWo^+(3s-)  +  sin2  0Wo^+{®n)  +  Wi,  (22) 

6  b 

P'2,2  =  WohA^^)  +  sin2  01^off+(®S-)  +  cos2  0Wo^+(®n)  +  Wi,  (22') 

Hz, 3  =  Woh:C'^~)  +  cos^tfWoirJ^E-)  +  OWohA^'H)  +  Wi,  (23) 

H4,4  =  WojjA^H)  +  sin^^Wo/rJ^S-)  +  cos^  0  WoH^("n)  +  Wi,  (23') 

Hz, 5  =  +  W2,  (24) 

where 

-  Ih.  (25) 

and 

W2  =  —  2  Ih  +  Ao,  (25') 

Using  the  basis  functions  (l9)-(2l)  one  obtains  the  following  non-zero  off-diagonal  ma¬ 
trix  elements: 

H,,,  =  sm9  cos«[lVo  (’")!,  (26) 

o  b 

Hz,a  =  sin^  cos0[Woji,(2e-)  -  WonA^H)],  (26') 
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1^3,5  =  sine  VoHt,  (27) 

=  -  cos  eVoHt,  (27') 

■^2.5  =  VoH^.  (27") 

2.3  Excited  state  (H20)+ 

The  (H20)'*‘  first  excited  state  configuration  is 

A^Ai  (Ia)^(2a)^(l6;,2)2(3a)^(l6y)^  (28) 

We  describe  the  electronic  structure  of  the  A  ^Ai  state  by  following  basis  functions: 

=  XixXiyXiyXizX2,  (29) 

^2  =  X\xXlyXiyXlzXz,  (29') 

^3  =  XixXiyXiyXizXiz^  (30) 

^4  =  XlxXlxXlyXiyXlz-  (30') 

The  functions  and  $2  correlate  with  the  asymptotic  0(®P)HaH^  and  0(^P)H+H6 


configurations,  respectively,  and  the  functions  $3  and  $4  correlate  with  the  asymptotic 


0  (2p)H+H^.  The  diagonal  matrix  elements  of  the  basis  functions  (29)-(30) 

are  as 

follows: 

^1,1  —  UoH„(^n)  +  WQ^+(®n)  -f  —  Iff, 

0 

(31) 

Hz, 2  =  ^oH+(^n)  +  —  Ijj, 

(310 

Hz,z 

=  ^o-H+(^n)  -b  sin^^^W'o-H+C^S)  +  cos2tfWo_^+(2n)  +  W2, 

(32) 

H4,4 

=  +  cos2  0Wo_^+(2s)  +  sin^eUo-jj^Cll)  +  W^. 

(32') 

Using  the  basis  functions  (29)-(30)  one  obtains  the  following  non-zero  olF-diagonal  ma- 

trix  elements: 

Hi,z  =  VoH^, 

(33) 

H2,Z  =  —  cos  e  VoHi,  5 

(33') 

-^2,4  =  sin^VoHj, 

(33") 

Hz, A  =  sintf  cos/?  [Wo_jj+(2s)  _  UQ.rj+[^'a)\ 

b 

(34) 
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In  the  absence  of  the  Ar  atom  the  X  and  A  states  of  (HjO)'*'  ionic  molecule  are 
uncoupled  with  one  another. 


2.4  (ArH20)‘*'  system 


As  it  has  been  assumed  by  us,  the  electronic  structure  of  Ar  and  Ar"*"  is  described 
in  terms  of  s-symmetry  X4  AO.  Within  the  spirit  of  the  DIIS  method  we  will  distin¬ 
guish  between  configurations  with  the  charge  located  on  the  H2O  molecule  and  with 
those  where  it  is  located  on  the  Ar  atom.  We  consider  first  the  Ar(H20)+  configuration 
states  which  correlate  asymptotically  with  the  configurations  ArOHaH^,  ArOH+Hj, 
and  ArO~H+H^.  The  basis  functions  (19)-(21)  and  (29)-(30)  present  all  possible  ori¬ 
entations  of  the  (H20)'*"  electronic  configuration.  Adding  the  AOs  of  the  Ar  atom  we 
obtain  nine  basis  functions  of  the  Ar(H20)‘'‘  configuration.  In  what  follows  these  basis 
functions  along  with  the  asymptotic  atomic  configurations  are  presented: 


ArOHaH+: 


ArOH+Hfc: 


ArO-H+H^ 


=  XliXlyXlj,XlzX2X4X4>  (35) 

^2  =  XlxXixXiyXizX2X4X4,  (35') 

^3  =  XlxXlyXl2XlaX2X4X4.  (35") 

^4  =  Xla:XlyXlyXl^X3X4X4)  (36) 

$5  =  XliXla:XlyXl2X3X4X4)  (36') 

^6  =  XlsXlyXl«XizX3X4X4.  (36") 

^7  =  XlxXlyXlyXl^Xl^X4X4.  (37) 

^8  =  Xla:XliXlyXl2XlzX4X45  (37') 

^9  =  XixXiiXiyXiyXi2X4X4-  (37") 


Next  the  Ar"^H20  configuration  and  the  basis  functions  which  correlate  with  the 
asymptotic  atomic  configurations  Ar+OHoHj,  Ar‘*‘0“HoH^,  and  Ar+0“H+H{,  are 
considered.  Since  an  arbitrary  3-dimensional  geometry  has  to  be  described  we  have 
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to  include  all  orientations  of  the  2p  AO  of  0  and  not  only  those  presented  by  the  H2O 
basis  functions  (3)- (5).  This  provides  (3  x  3  =)9  basis  functions  of  Ar"^H20.  We  have 
to  distinguish  also  between  the  states  with  different  atomic  spin  orientations  in  H2O 
relative  to  the  Ar"^  spin  which  is  assumed  to  be  5^  =  -.  As  a  result,  the  number 
of  diabatic  electronic  configurations  is  doubled  so  that  the  number  of  basis  functions 
becomes  (9x2  =)18.  We  shall  present  below  some  of  the  Ar“^H20  basis  functions 
omitting  those  of  them  with  expressions  that  can  be  easily  derived.  The  configurations 
with  positive  and  negative  0  atom  spin  orientation  are  described  by  with  odd  and 

even  k,  respectively. 

Ar+OHaHfc 

=  XlxXlt/Xl!,XlzX2X3X4, 

(38) 

^2  =  XliXlyXl!/XlzX2X3X4, 

(380 

^3  =  XlxXliXl!/XlzX2X3X4, 

(38") 

$5  =  XlxXl!/XlxXlzX2X3X4, 

(38"') 

Ar+O-HaH^ 

^7  =  XlxXlyXlyXlxXlzX2X4, 

(39) 

^8  =  XlxXlyXlyXlxXlzX2X4, 

(390 

^9  =  XlxXizXlyXlaXizX2X4, 

(39") 

$11  =  XlxXlxXlyXiyXlxX2X4, 

(39'") 

Ar+O-H+Hft 

$13  =  XlxXlyXlyXlxXlzX3X4, 

(40) 

$18  =  XlxXlxXlyXls,XlzX3X4, 

(40') 

The  ofF-diagonal  matrix  elements  between  the  states  with  different  oxygen  AO 
orientation  coincide  with  those  (see  Eqs.(15)-(17)  and  (26)-(27))  of  the  (H20)+  and 
H2O  molecules.  The  coupling  between  the  Ar‘*‘H20  states  with  different  AO  spin 
orientation  (between  and  $2  of  Eqs.(38),  for  example)  is  taken  into  account  by  the 
states  of  the  0~H  fragment.  We  will  treat  this  coupling  in  a  way  similar  to  the  way 
the  H'^H  —  HH'*'  coupling  was  treated  (see  Sec.2.1).  The  0~H  adiabatic  potentials 
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are  presented  as  the  eigenvalues  of  Eq.(lO).  In  the  0“H  case  the  diabatic  states  are 
distinguished  by  the  atomic  spin  orientation 


=  XiwXiwXiK>Xi^'Xi<TX2y  (41) 

h  =  XlKXl^Xlir>Xi^>XlaX2-  (41') 


The  coupling  between  the  diabatic  electronic  configurations  of  Ar+H20  (Eqs.(38) 
-(40))  and  Ar(H20)+  (Eqs.(35)-(37))  is  tahen  into  account  by  the  state  of  the 

(ArH)+  fragment.  The  adiabatic  potentials  of  the  (ArH)+  fragment  are  the  eigenvalues 
of  the  equation 


UArH+  -  E  Vf^ArH)+ 
^(>lrH)+  Uat+H-E 


(42) 


and  consequently  17^;.^+  and  U^r+H  are: 


UArH+  —  -  [I^(ArH)+ (-X”  ^S)  -f-  W'’(ArH)+ (2  ^S)  —  Wd],  (43) 

Uat+H  =  ^\W(ArH)+{X^i:)  +  W^(ArH)+ (2  ^S)  +  W^],  (43') 

where 

=  ^[W'(ArH)+(2lS)  -  W^(^,h)+(X1E)]2  -  (44) 


3.  Diatomic  potentials 

The  required  adiabatic  potentials  to  evaluate  matrix  elements  can  be  taken  from 
empirical  studies  or  ab  initio  calculations.  In  our  treatment,  however,  most  of  the 
excited  diatomic  potentials  are  considered  as  fitting  functions  with  adjustable  parame¬ 
ters.  While  doing  that  we  keep  as  much  as  possible  their  character  and  generally  known 
features  from  empirical  studies  or  ab  initio  calculations. 

The  diatomic  potentials  of  the  ionic  molecules  of  the  type  A'^B  or  A~’B^  are 
presented  as  a  sum  of  a  short  range  valence  potential  and  a  long  range  electrostatic 
potential  U^i 

U  =  +  U,u  (45) 

In  the  case  of  an  asymptotic  configuration  with  one  charged  atom,  for  example 
the  electrostatic  interaction  is  of  a  polarization  origin 

Uei  =  P{R)  =  -OB  7(-R)/i2‘‘,  (46) 


10 


where  R  is  the  interatomic  separation,  ag  is  the  polarization  coefficient  of  the  neutral 
atom  B,  and  7  is  a  dumping  term  introduced  to  prevent  the  unphysical  behaviour  of 
P(R)  at  small  separations  R.  The  suggested  form  for  the  functional  dependence  7(R) 
is® 

'^{R)  =  Ry[R^^  +  {rA+ +  (47) 

with  r^+  and  tb  being  the  atomic  radii.  A  suggested  form  for  the  attractive  A^B~ 
Coulomb  potential  is 


U,i  ^  U,  =  -Uo/[R^  +  irA^  (48) 

where  Uq  is  equal  to  14.4eVA.  The  term  (rx+  +rB-)'*  in  Eq.(48)  prevents  the  unphysical 
behaviour  of  Uc  at  small  separations  R. 

3.1  O-H  potentials 

The  matrix  elements  (6)-(9)  of  H2O  and  (22)-(24)  and  (3l)-(32)  of  (H20)+  include, 
in  particular,  adiabatic  and  diabatic  potentials  of  the  diatomic  fragments  OH,  (OH)"'', 
and  {OH)“.  Since  the  diabatic  potentials  are  expressed  in  terms  of  the  adiabatic  ones 
(see  Eqs.(ll),(l3))  we  will  discuss  the  adiabatic  potentials  only. 

The  OH  potentials  are:  X^H,  3  ^H,  C  1  The  OH  ground  and  excited 
states,  i.e.  X^H  and  1^E“  respectively  have  the  atomic  ground  states  0(®P)  and  H(^<S') 
as  dissociation  limit.  The  OH  excited  states  3  ^H  and  have  0~{^P)  and  H"*"  as 
dissociation  limit.  Except  for  the  state  which  is  repulsive  the  other  three  terms  are 
bound. The  ground  state  X^H  potential  is  described  in  terms  of  a  modified  Morse 
potential^®  which  fits  the  ab  initio  results  of  Ref.23.  The  OH  excited  state  potentials 
are  considered  as  fitting  functions  with  adjustable  parameters.  The  potentials  3  ^H 
and  which  correlate  asymptotically  with  0“+H+  include  the  attractive  Coulomb 

term  (48).  As  it  was  mentioned  above  (see  Eq.(l2))  that  the  ^H  diabatic  potentials  are 
coupled  with  each  other  and  therefore  an  exchange  term  Vn  is  required  which  will  be 
considered  as  a  fitting  function. 

The  (OH)"*"  adiabatic  potentials  are:  X®E“  and  A^H.  Both  are  bound.  The  ground 
state  X  potential  which  has  the  0(®P)  +  H"^  as  dissociation  limit  is  presented  as 
a  sum  of  a  Morse  potential  and  the  polarization  term  (46)  (it  fits  reasonable  well  the 
ab  initio  results  of  Ref.23).  The  excited  state  A  potential  has  0'''(^5)  +  H  as 
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dissociation  limit^^  and  is  also  presented  as  a  sum  of  a  Morse  potential  (with  adjustable 
parameters)  and  the  polarization  term  (46). 

The  (OH)  potentials  are:  X  2  ^E*^  and  ^H.  All  of  these  states  have  the 
O  (^P)+H  as  the  dissociation  limit.  The  (OH)~  ground  state  which  is  bound^^  is 
presented  as  a  sum  of  a  Morse  potential  and  a  polarization  term  (46).  The  dissociation 
energy  of  the  X  state  is  considered  as  an  adjustable  parameter.  In  contrast  to  the 
bound  state  the  excited  states  are  repulsive. 

3.2  H-H  potentials 

The  H2O  and  (H20)+  matrix  elements  include  the  excited  ^E^  H2  potential,  the 
ground  ^Ey  and  excited  ^E^  potentials  of  and  the  ^E  potential  of  the 
repulsion.  The  repulsive  ^E^  potential  is  presented  in  terms  of  a  modified  anti-Morse 
potential^o  which  fits  the  ab  initio  results.^^  The  ^E^  and  ^Eu  potentials  of  HJ  are 
assumed  to  be  Morse  and  anti-Morse  potentials,  respectively,  which  fit  the  results  of  aji 
accurate  calculation.^® 


3.3  Ar-H  potentials 

The  (ArH20)“^ matrix  elements  include  a  van  der  Waals  (vdW)  potential  X^E  of 
ArH,3®  three  (ArH)+  potentials  of  the  type  X  ^E,  2  ^E,  and  ^E,  and  a  1  Coulomb 
potential  of  Ar'^H"^.  In  the  (ArH)^  molecule  the  ground  state  X  ^E  has  the  Ar+H^ 
as  the  dissociation  limit  whereas  the  excited  states  2  ^E  and  ^E  dissociate  into  Ar++H. 
Since  we  replaced  the  Ar"^  3p-AO  by  an  s-AO  (see  Ch.2)  we  cannot  distinguish  be¬ 
tween  the  (ArH)+  and  Ar+H+  potentials  and  the  known  adiabatic  potentials  of  these 
molecules.^^-^^  We  assume,  however,  that  the  approximated  ground  state  X^E  poten¬ 
tial  coincides  with  the  real  (ArH)*^  ground  state  potential.  This  potential,  presented 
as  a  sum  of  a  Morse  potential  and  the  Ar+H  polarization  term  (46),  fits  the  known  ab 
inxtio  potential.^®’^^  The  excited  state  potentials  of  (ArH)"*”  are  presented  as  a  sum  of 
an  anti-Morse  potential  and  the  H+Ar  polarization  term  (46).  The  Ar+H+  potential  is 
presented  in  terms  of  a  repulsive  Coulomb  potential. 

3.4  Ar-O  potentials 

The  matrix  elements  include  the  following  Ar-0  diabatic  potentials:  the  ^E”  and 
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of  ArO,  the  '‘S",  ^n,  and  of  Ar+0,  the  2e+  and  of  ArQ-,  and  the  ®n, 
^n,  and  of  Ar'^O".  Unfortunately,  no  ab  initio  calculations  of  Ar-O  interac¬ 
tions  are  known  to  us.  Some  indirect  information  can  be  obtained  from  the  empirical 
study  of  ArO  diatomic  in  the  excited  ^  E"^  state^^  and  from  semiempirical  calculation  for 
the  ArOH  system.®®  A  general  idea  about  the  behaviour  of  the  argon-oxygen  potentials 
in  different  states  can  be  obtained  from  ab  initio  calculations  of  diatomics  with  similar 
electronic  structure.  For  instance,  the  ArO~  electronic  structure  is  similar  to  that  of 
the  well  studied  ArF  molecule®®  and  the  Ar'*'0“  electronic  structure  is  similar  to  that 
of  FCl.®^ 

Due  to  lack  of  Ar-0  potentials  (and  their  ions)  these  had  to  be  constructed  .  In 
doing  that  we  took  into  accoimt  the  Ar  and  O  radii,  as  well  as  some  indirect  informa¬ 
tion  mentioned  before.  The  ion-atom  and  ion-ion  potentials  include  the  electrostatic 
terms  (46)  and  (48),  respectively.  It  is  important  to  note  that  axgon-oxygen  potentials 
parameters  cannot  be  treated  as  adjustable  parameters,  like  it  was  done  in  the  case  of 
some  of  the  0-H  potentials,  since  there  is  not  available  any  data  to  be  employed. 


4.  3-center  term 

In  order  to  improve  the  results  (see  the  following  Section)  we  incorporated  into 
the  (H20)"*"  matrix  elements  the  so  called  3-center  term.  26,38  (H20)"*’  case  the 

suggested  3-center  term  is 


^3c  =  —gexp[—S]{l  +  bsin^6)  (49) 

where 

S'  =  a  [(iii.2  -  Ro)^  +  {Ri,3  -  Ro)%  (49') 

0  is  the  HOH  angle,  Ri^2  and  Ri^s  are  the  0-H  distances.  The  parameters  of  Eqs.(49) 
are  as  follows:  g  —  1.2eV,  b  =  0.22,  a  =  3A~®,  Rq  =  1.03A. 


13 


5.  The  results  of  H2O  and  (H20)"’‘  calculation 

The  adjustable  parameters  of  the  excited  diatomic  potentials  (Ch.3)  and  the  3- 
center  term  (Ch.4)  were  used  to  fit,  at  leaist  approximately,  the  known  empirical  and/or 
ab  initio  parameters  of  the  H2O  molecule  and  the  (H20)'’’  ion. 

The  results  for  the  ground  state  X  ^Ai  H2O  molecule  calculation  are  presented 
in  Table  1.  These  are  compared  with  the  H2O  empirical  parameters  as  well  as  with 
the  parameters  obtained  from  the  DIM  and  the  ab  initio  calculations.  It  is  noticed  that 
our  calculation  provides  an  H2O  geometry  which  is  very  close  to  the  empirical  one. 
In  addition  our  vibrational  frequencies  are  in  reasonable  agreement  with  the  empirical 
and  the  ab  initio  frequencies.  On  the  other  hand  o\ir  calculated  dipole  moment  is  about 
25%  lower  than  the  empirical  value  so  that  the  Ar”^— H2O  electrostatic  interaction  is 
underestimated.  Also,  our  calculated  H2O  energy  is  higher  than  the  empirical  value 
by  about  0.66eV.  This  may  look  like  a  significant  discrepancy  however  we  are  interested 
in  the  energy  gaps  between  the  different  states  and  not  in  the  absolute  values.  These 
energy  gaps,  as  it  will  be  shown  below,  are  of  much  higher  accuracy. 

The  results  of  the  ground  state  X'^Bi  (H20)+  calculation  are  presented  in  Table 
2.  The  calculated  (H2  0)'*"  energy  is  higher  than  the  empirical  value  by  0.56eV  but 
the  difference  between  the  (H20)+  and  H2O  energies  is  reasonable  close  to  the  H2O 
ionization  potential. 

As  was  mentioned  in  the  Introduction  the  (H2  0)+  is  formed  in  the  excited  elec¬ 
tronic  state  A  ^Ai.  The  various  calculated  parameters  related  to  the  A  ^Ai  (H20)+ 
are  presented  in  Table  3.  Like  in  the  ab  initio  calculations'*®’^®  the  equilibrium  of  the 
(H20)+  geometry  in  the  A  state  is  collinear  (in  this  geometry  the  A  PES  is  tangent  to 
the  ground  state  of  the  X  PES)  but  the  equilibrium  bond  length  is  found  to  be  some¬ 
what  larger  than  the  ab  initio  value.  The  following  parameters:  the  X  — >  A  excitation 
energy,  the  excess  energy  AE  related  to  the  X  A  transition  (for  the  A  equilibrium 
energy),  and  the  vibrational  frequencies  w,  all  lie  in  a  reasonable  proximity  of  the  cor¬ 
responding  parameters  of  the  ab  initio  calculations^®’"*®  (as  well  as  with  these  of  the 
empirical  model"*®).  The  exceptional  case  is  the  bending  frequency  which  fits  the  one 
by  Brommer  at  al."*®  but  differ  significantly  from  the  one  by  Smith  at  al."*® 

The  A  state  (H20)+  ion  formed  in  the  charge  transfer  reaction  (l)  is  expected 
to  be  found  in  highly  excited  vibrational  bending  states  (0,V2,0).^  Since  the  vibrational 
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bending  levels  play  an  important  role  in  the  study  of  the  charge  transfer  reaction  we 
calculated  them  for  the  A  PES.  These  calculations  were  done  for  the  overall  rotational 
state  K=0.  The  calculated  vibrational  modes  are  presented  in  Table  4. 

In  order  to  study  the  stretching-bending  coupling  effect  on  the  vibrational  bending 
modes  calculations  were  done  twice.  In  the  first  (column  I  in  Table  4)  the  stretching 
motion  is  restricted  to  the  zero  mode  vx  =  U3  =  0  so  that  the  stretching— bending 
coupling  effect  is  of  a  negligible  importance.  According  to  this  calculation  the  energy 
gaps  hu;  between  the  bending  levels  are  almost  constant  for  low  bending  levels  but 
increase  once  V2  >  4,  In  other  treatments  of  the  bending  modes'^®”'^®  the  energy  gaps 
increase  with  V2  from  the  very  begiiming. 

In  the  second  treatment  (column  II  in  Table  4)  the  stretching  modes  as  well  as 
the  stretching-bending  couplings  were  taken  into  account.  While  comparing  the  results 
of  the  two  calculations  (Table  4)  it  is  noticed  that  the  main  effect  of  the  coupling  is 
in  lowering  the  bending  levels.  This  effect  is  stronger  the  larger  is  V2.  For  instance, 
the  (0,1,0)  level  is  almost  the  same  according  the  two  calculations  whereas  the  (0,9,0) 
level  becomes  lower  by  0.22eV  according  the  second  calculation.  Due  to  the  levels  shift 
caused  by  the  stretching-bending  coupling  the  energy  gaps  between  the  various  bending 
states  remain  almost  constant  along  the  whole  studied  interval,  i.e.  2  <  V2  <  12.  Thus 
according  to  our  treatment  the  stretching-bending  couplings  affect  significantly  the 
bending  levels  manifold  and  it  follows  that  the  bending  levels  obtained  within  the  rigid 
OH  bond  approximation  are  far  from  being  accurate. 

6.  (ArH20)+  potential  energy  surfaces 

A  charge  transfer  reaction  can  be  treated  within  the  electrostatic  approximation  if 
the  electron  transition  occurs  at  intermolecular  distances  exceeding  the  strong  coupling 
range.  Using  the  PESs  of  the  present  calculation  it  is  possible  to  estimate  this  range 
and  consequently  to  determine  the  limits  of  the  electrostatic  approximation. 

Let  us  consider  first  the  dependence  of  the  (ArH20)+  PESs  on  the  HOH  valence 
angle  1?  (Fig.2).  These  PESs  were  obtained  for  Ar+  ion  approaching  one  of  the  H  atoms 
along  the  0-H  line  (z-axes  of  Fig.l).  The  OH  distances  of  the  H2O  molecule  are  fixed 
and  equal  to  the  OH  equilibrium  distance  (0.958A).  The  PESs  of  Fig.2  demonstrate  a 
typical  avoided  crossing  between  two  electronic  states  at  «  108°  - 110°,  close  to  the 
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equilibrium  angle  of  104*8°.  Since  the  region  of  the  avoided  crossing  is  accessible  in 
the  vibrational  ground  state  the  bending  motion  may  induce  the  charge  transfer.  The 
avoided  crossing  gap  AE  increases  when  the  oxygen-argon  distance  RoAr  decreases  and 
becomes  comparable  to  the  electrostatic  potential  at  the  distance  Rqat  ^  SA. 

The  adiabatic  PESs  dependence  on  the  argon  position  along  a  straight-line  trajec¬ 
tory  (see  Fig.3)  is  shown  in  Fig.4afor  a  fixed  (equilibrium)  H2O  geometry.  The  PESs  are 
presented  for  three  electronic  states,  namely  k=2,4,5.  (The  ground  state  k=l  is  omitted 
as  it  asymptotically  correlates  with  (H20)'^  in  the  ground  state  X  and  the  state  k=3 
is  ignored  as  it  is  not  coupled  with  the  states  of  interest  due  to  a  different  symmetry.) 
The  charge  delocalization  caused  by  the  electronic  coupling  is  demonstrated  in  Fig ,4b. 
According  to  the  results  of  the  calculation  the  effect  of  the  coupling  is  negligible  small 
at  distances  |p|  >  SA  but  becomes  relatively  strong  when  |p|  «  3A  which  leads  to  the 
deformation  of  the  potential  curves  (Fig.4a)  and  to  charge  delocalization  (Fig.4b).  Due 
to  this  coupling  the  k=2  potential  curve  of  Fig.4a  (which  correlates  cusymptotically  with 
Ar^+H20)  decreases  as  \p\  decreases  and  the  k=4  potential  (which  correlates  asymp¬ 
totically  with  Ar+(H20)“*"  )  increases.  Since  the  charge  delocalization  takes  place  when 
the  states  k=2,4  start  to  diverge  the  charge  is  located  on  both,  on  Ar  and  on  H2O  . 
Whereas  around  p=0  the  k=5  electronic  state  becomes  the  one  with  the  largest  Ar 
charge  the  opposite  is  encountered  for  the  k=4  state  where  the  charge  is  mainly  located 
on  H2O  .  The  regions  where  the  charge  distribution  is  most  affected  (those  are  the 
regions  of  the  avoided  crossing)  are  those  in  the  vicinity  of  \p\  «  3.5A  (2  — 4)  and 
|p|«iA(4^5). 

For  the  results  presented  in  Fig.4  (as  well  as  those  in  Fig.2),  the  range  of  the 
strong— coupled  region  was  found  to  be  roughly  equal  to  4A.  Comparing  this  value  with 
the  reactive  interaction  range  Rj-e  it  is  possible  to  conclude  whether  the  electrostatic 
approximation  is  valid  or  not.  We  suggest  to  determine  Rre  as  the  average  impact 
parameter 


where  a  is  the  reactive  cross  section.  At  near-thermal  energies  the  cross  section  of  the 
Eq.(l)  charge  transfer  reaction  is  very  large,  i.e.  SO-fiOA^.^’'  The  respective  reactive 
interaction  range  of  Eq.(50)  is  about  4 A,  the  same  as  the  strong-coupling  range.  The 
conclusion  can  be  made  that  the  near-thermal  reaction  can  be  treated  most  probably 
by  the  electrostatic  approximation.  At  the  suprathermal  collision  energies  of  1-lOeV 
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the  reaction  cross  section  is  12-18 The  respective  Rre  values  are  2-2.4 A.  It  follows 
that  in  this  energy  interval  the  electron  hop  occurs  mainly  in  the  strong— coupling  region 
so  that  the  description  of  the  reaction  in  the  terms  of  electrostatic  interactions  is  not 
always  sufficient. 


7.  Charge-transfer  reaction  study 

The  charge  transfer  reaction  is  expected  to  be  governed  by  the  complex  formation 
at  low  collision  energies  and  by  a  direct  mechanism  at  suprathermal  energies.'^®  In  the 
forthcoming  theoretical  treatment  we  consider  the  suprathermal  collision  energy  region. 
Also  in  what  follows  we  will  refer  to  the  (H20)"^  in  the  first  excited  electronic  A  state 
exclusively. 

The  H2O  energy  in  the  vibrational  ground  state  (0,0,0)  is  -8.83eV  (Table  1)  and 
the  Ar  ionization  potential  is  15.76eV  so  that  the  energy  of  the  reagents  in  reaction  (l)  is 
6.93eV.  The  energy  of  the  products  (the  Ar  and  the  A  (H20)+  in  the  vibrational  grormd 
state)  is  4.99eV  (Table  4),  so  that  the  energy  difference  is  equal  to  1.94eV.  This  excess 
energy  is  converted  into  the  internal  energy  of  the  (H20)'^  ion.^’^^  However,  as  the 
equilibrium  geometries  of  the  reagent  H2O  and  the  product  (H20)+  differ  significantly 
in  their  valence  angles  (104.8°  versus  180°)  the  bending  (and  not  the  stretching)  modes 
are  the  ones  expected  to  be  excited.  This  expectation  was  confirmed  experimentally.^ 

In  the  energy  interval  1  <  £■  <  lOeV  the  collision  time  is  of  the  order  of  4  •  10“^^  - 
2  ■  10  ^®sec.  The  duration  of  the  H2O  rotation  is  usually  longer  and  we  may  assume 
the  fixed  rotor  approximation.  The  collision  time,  however,  is  much  longer  than  the 
H2O  bending  mode  period  which  is  ~  2  •  10~^^sec.^®  Consequently  the  H2O  molecule 
has  the  time  to  cross  more  then  once  the  avoided  crossing  region  (Fig. 2)  with  some 
probability  to  form  the  (H20)+  in  the  vibrationally  excited  (0,t;2,0)  state.  This  state 
can  be  considered  as  being  near-resonant  to  the  ground  state  level  of  the  reactants 
Ar"'‘-1-H20  system. 

At  large  distances  the  reactants  Ar+-|-H20  ground  state  (E=6.93eV)  overlaps  with 
the  (0,9,0)  vibrational  bending  state  of  the  products  Ar+-f-(H20)+  (E=6.93eV,  see 

Table  4  column  11) .  In  the  absence  of  the  stretching-bending  coupling  (Table  4  column 
I)  the  reactants  ground  state  is  in  near  resonance  with  the  products  (0,8,0)  level,  in 
accordance  with  the  model  calculations  of  Ref.l.  In  the  interaction  region  of  the  two 
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molecular  species  the  energies  of  both,  reagents  and  products,  are  shifted  so  that  the 
(0,0,0)  Ar“*'H20  vibrational  level  may  overlap  with  the  Ar(H20)'^  (0,V2,0)  levels  whose 
V2  is  different  from  the  one  found  in  the  cisymptotic  region.  In  order  to  study  the 
behaviour  of  the  Ar'^H20  and  Ar(H20)'^  levels  we  performed  two  kinds  of  calculations. 
In  both  we  assumed  that  the  argon  ion  (atom)  is  approaching  (and  receding)  the  water 
molecule  along  a  straight-line  trajectory. 

7.1  The  simplified  treatment 

In  the  simplified  treatment  it  is  assumed,  like  in  Ref.l,  that  the  geometry  of  the 
water  molecule  (ion  )  is  not  affected  by  the  Ar  ion  (atom)  and  that  the  system  energy 
is  equal  to  the  sum  of  a  PES  W  and  a  vibrational  energy  Ey  of  an  isolated  molecule 

E{p',Ki,V2)  =  W{p-,b,^)  +  E^{v2)  (51) 

In  Eq.(5l)  /s  is  a  coordinate  along  the  Ar  trajectory  (Fig. 3),  b  is  the  impact  parameter, 
and  7  comprises  the  intermolecular  orientation.  The  diabatic  energy  levels  along  a 
straight-line  trajectory  at  an  impact  parameter  b=3A  are  shown  in  Fig. 5  for  three 
different  H2O  orientations.  Fig.5  presents  energy  levels  of  both,  reagents  Ar‘*'H20  and 
products  Ar(H20)+  systems.  In  the  case  of  the  Ar+H20  system  W  is  the  PES  of  the 
frozen  (equilibrium)  H2O  and  is  the  vibrational  energy  of  the  isolated  H2O  and  in 
the  case  of  Ar(H20)+  system  W  is  the  PES  of  the  frozen  (equilibrium)  (H20)+  and 
E„  is  the  vibrational  energy  of  isolated  (H20)+.  The  (H20)+  Ey  results  follow  for  the 
exact  calculation  (see  Table  4  column  II). 

It  is  seen  in  Fig.5  that  at  large  distances  the  Ar(H20)+  levels  form  a  flat  manifold 
for  all  (H20)+  orientations.  When  the  distance  becomes  relatively  small  the  levels 
are  slightly  decreasing  due  to  the  attractive  ion-induced  polarization  forces  which  are 
isotropic.  The  behaviour  of  the  reagents  Ar+H20  level  is  different  due  to  charge- 
dipole  forces  which  are  much  stronger  than  the  polarization  forces  and  can  be  either 
attractive  or  repulsive  depending  on  the  H2O  orientation.  In  the  case  of  Fig.5a  (the 
relevant  orientation  is  shown  in  Fig.3)  the  charge-dipole  interaction  is  repulsive  and 
consequently  the  Ar+H20  level  is  increasing  when  |/»|  decreases.  In  the  vicinity  of 
p  =  0  the  direct  repulsion  between  the  Ar"*"  and  H  contributes  to  a  further  increase  of 
the  Ar+H20  energy.  As  a  result  the  reagents  Ar''"H20  level  crosses  four  Ar(H20)+ 
levels.  In  case  of  Fig.5b  for  the  presented  orientation  the  charge-dipole  interaction 
is  attractive  but  in  the  vicinity  of  p  =  0  this  electrostatic  interaction  is  compensated 
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by  the  direct  interatomic  repulsion.  There  are  no  level  crossings  in  this  case.  Fig. 5c 
presents  the  case  of  a  different  interactions  for  incoming  and  outgoing  trajectories.  In 
particular,  the  charge-dipole  interaction  is  attractive  for  />  <  0  and  repulsive  for  p  >  0. 
In  the  vicinity  of  p  =  0,  because  of  the  proximity  of  the  argon  atom  to  one  of  the 
hydrogen  atoms,  a  strong  repulsion  is  demonstrated.  The  Ar'^H20  ground  vibrational 
level  of  Fig. 5c  crosses  a  few  times  the  Ar(H20)"^  excited  levels.  According  to  Fig. 5 
the  products  (1120)“^  which  is  formed  along  the  b— sA  trajectory  is  formed  in  different 
vibrational  bending  (0,t;2,0)  states,  namely  V2  =  8,9,10,11.  We  do  not  consider  here 
the  trajectories  with  the  impact  parameters  b  <  sA  as  the  straight-line  assumption  is 
most  probably  not  valid  anymore. 

7.2  Full  scale  calculation 

In  contrast  to  the  simplified  calculation,  in  the  full  scale  calculation  the  H2O  and 
(H20)'^  vibrational  modes  are  calculated  not  for  isolated  species  but  in  the  field  of 
the  argon  ion  (or  atom).  In  Fig. 7a  are  presented  the  vibrational  eigenstates  of  the 
Ar'*"H20  and  the  Ar(H20)'^  systems  as  a  function  of  p  (see  Fig.6)  and  fixed  value 
of  b  (=3 A),  calculated  employing  average  PESs  obtained  following  a  rotation  around 
the  x-axis  (see  Fig.6).  The  two  PESs  used  in  the  calculation  were  constructed  for  the 
two  relevant  adiabatic  surfaces  (obtained  by  the  diagonalization  of  the  DIIS  matrix) 
which  then  were  chosen,  for  the  sake  of  this  calculation,  to  be  diabatic.  These  two 
PESs  become,  once  p  ^  00,  the  Ar~*“-|-H20  and  Ar+(H20)'^  surfaces.  It  was  found 
that  this  procedure  holds  as  long  p  >  2 A,  The  effect  of  the  argon— water  interaction 
on  the  charge  distribution  and  HOH  equilibrium  configuration  is  shown  in  Fig.Tb  for 
the  electronic  state  correlated  asymptotically  with  Ar+H20.  While  the  Ar+  approaches 
H2O  the  charge  concentrates  on  the  HOH  fragment  so  that  at  relatively  small  distances 
(p  <  4 A)  the  electronic  configuration  of  the  system  is  rather  described  as  Ar(H20)‘^ 
instead  of  Ar‘^H20.  The  increase  of  the  HOH  charge  is  also  accompanied  by  the  increase 
of  the  HOH  valence  angle. 

A  different  situation  is  encountered  in  case  of  the  asymptotic  Ar(H20)‘^  system. 
Here  the  Ar  atom  remains  uncharged  up  to  small  values  of  p.  Only  in  the  vicinity  of 
p  =  0  the  Ar  atom  acquires  small  quantities  of  charge  (about  20%).  The  approaching 
Ar  atom  does  not  affect  on  the  (H20)‘^  equilibrium  geometry  which  remains  collinear 
throughout  its  approach. 
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The  effect  of  coupling  on  the  (ArH20)'’'  energy  levels  is  shown  in  Fig.Ta.  At  large 
distances  [p  >  SA)  where  the  electrostatic  forces  are  dominant  the  Ar(H20)"'"  levels 
(0,V2,0)  form  a  flat  manifold,  like  in  the  simplified  treatment  (see  Fig.5).  The  same 
situation  is  encountered  in  the  case  of  Ar"*'H20  (0,0,0)  system  (this  is  in  contrast  to 
what  was  observed  in  the  simplified  model  —  see  Fig.5).  At  small  p  values,  in  the  strong- 
interaction  region,  all  levels  decrease.  It  is  important  to  emphasize  that  the  Ar(H20)'*' 
(0,9,0)  level  is  in  near  resonance  with  the  Ar'*‘H20  ground  state.  In  this  sense  our  more 
extensive  treatment  differs  significantly  from  the  simplified  calculation. 

Similar  results  were  obtained  for  other  orientations  of  the  HOH  rotation  axis  and 
therefore  they  will  be  not  presented  here. 


8.  Conclusions 

1.  The  semiempirical  diatomics-in-ionic-systems  (DIIS)  method  proves  to  be  an 
effective  tool  for  treating  the  PESs  of  the  (ArH20)+  system.  Within  the  DIIS  method, 
like  in  the  DIM,  the  Hamiltonian  matrix  elements  are  expressed  in  terms  of  diatomic 
potentials  but  in  the  DIIS  along  with  adiabatic  diatomic  potentials  one  employs  also 
the  diabatic  potentials  which  describe  electronic  configurations  with  localized  charges. 

2.  According  to  this  treatment  the  product  (H20)+  ion  (see  Eq.(l))  formed  in 
the  A  ^Ai  state  is  highly  vibrationally  (bending)  excited.  Also  it  was  found  that  the 
calculated  vibrational  bending  manifold  is  strongly  affected  by  the  stretching-bending 
coupling. 

3.  It  is  shown  that  the  electrostatic  approximation  may  be  applied  to  the  near- 
thermal  reactions  only.  At  suprathermal  collision  energies  this  approximation  is  not 
valid  anymore  and  the  electron  hop  occurs  mainly  in  the  strong  coupling  region.  Also 
it  is  established  that  the  electron  transfer  is  induced  by  the  H2O  bending  motion  and 
it  takes  place  at  the  region  of  Ar+H20 — Ar(H20)+  avoided  crossing. 

4.  The  Ar'’'H20  and  the  Ar(H20)+  vibronic  bending  modes  were  calculated  by 
two  different  ways,  a  simplified  version  and  a  full  scale  one.  In  the  simplified  version 
the  charge  is  assumed  to  be  located  either  on  the  Ar  or  the  H2O  and  the  effect  of  Ar+ 
(or  Ar)  on  H2O  (or  (H20)+)  is  ignored.  Consequently  it  was  found  that  the  (0,0,0) 
state  is  in  near  resonance  with  the  (0,9,0)  state  in  the  asymptotic  region  only  but  then 
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crosses  the  Ar+(H2  0)'^  manifold  as  the  distance  decreases.  The  two  assumptions  were 
not  included  in  the  full  scale  versions  and  as  a  result  we  found  that  (0,0,0)  state  of 
Ar^+HsO  is  in  near  resonance  with  the  (0,9,0)  bending  state  of  the  Ar+(H20)‘^ 
along  the  whole  interaction  region. 
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Table  1.  H2O  properties 


present 

calc. 

Empirical 

ab  initio®)  ab  initio’’)  DIM^  DIM®) 

OH  bond  length  (A) 

0.958 

0.957^) 

0.943  0.982  0.999  0.979 

HOH  angle  (deg.) 

104.8 

104.5^) 

106.7  102.3  10.9  100.3 

Dipole  moment  (D ) 

1.36 

1.87^) 

2.21 

Energy  (eV) 

-9.412 

-10.069’’) 

-8.812  -9.868  -9.37 

Energy  Eq  (eV) 

-8.831 

-9.605^) 

(010) 

1832 

1648®)  1598^’) 

1751 

Vibr.  CO  (cm-b(OOl) 

3518 

3832®)  3656^) 

4164 

(100) 

3898 

3943®) 

4261 

%ef.  39  %ef.  40  41  %ef.  42  ®)Ref.  43  %ef.  22  S^Ref.  21 


Table  2.  (H20)^  ground  state  properties 


Present 

calc. 

Empir.  ab  initio^^ 

OH  bond  length  (A) 

1.028 

(Ot) 

1.001  '  1.000 

HOH  angle  (deg.) 

112.2 

108.9^)  109.1 

Energy  E^  (eV) 

3.011 

H2O  ioniz.  pot.  I  (eV) 

12.406 

12.62^> 

Energy  Eg  (eV) 

3.575 

3.015^') 

(010) 

1711 

1388^^  1477 

Vibr.  0)  (era"*)  (001) 

3543 

3220^^  3381 

(100) 

3648 

3454^)  3438 

Ref.  44  Ref.  39  Calculated  as  Eq  + 1  of  H2O  (for  Eg  see  Table  1). 


Table  3.  (H20)‘'’  excited  state  A  properties 


Geometry 

ab  initio^^ 

ab  initio^^ 

Empir. 

model^^^ 

Equilibrium 

OH  bond  length 

(A) 

1.008 

0.978 

0.989 

(H20)+  a 

Energy  E^ 

(eV) 

4.370 

(010) 

1880 

799 

1501 

1679 

Vibr.  0)  (cm'^) 

(001) 

3640 

3648 

(100) 

4156 

3825 

Equilibrium 

Energy  E 

(eV) 

AE  =  E  -  E, 
e 

(eV) 

1.150 

1.61 

1.076 

(H20)+X 

X  — >  A  excit.  en. 

(eV) 

2.510 

2.59 

2.214 

Equil .  H2O 

Energy  E 

(eV) 

6.712 

%ef.  43  %ef.  45  ‘^^Ref.  46 


Table  4.  Vibronic  bending  levels  (0,V2,0)  of  (H20)'''  in  the  A  state  (in  eV). 


^1  ^^2 

la) 

energy 

levels  hco 

Ilb) 

energy 

levels  hco 

hco^) 

hco'^) 

0  0  0 

0  1  0 

4.995 

5.235  0.240 

4.993 

5.226  0.233 

0.208 

0.186 

0  2  0 

5.474 

0.239 

5.451 

0.225 

0.224 

0.207 

0  3  0 

5.712 

0.238 

5.672 

0.221 

0.235 

0.225 

0  4  0 

5.952 

0.240 

5.890 

0.218 

0.243 

0.234 

0  5  0 

6.194 

0.242 

6.107 

0.217 

0.249 

0.243 

0  6  0 

6.441 

0.247 

6.326 

0.219 

0.254 

0.249 

0  7  0 

6.691 

0.250 

6.545 

0.219 

0.258 

0.255 

0  8  0 

6.947 

0.256 

6.765 

0.220 

0.261 

0.258 

0  9  0 

7.205 

0.258 

6.984 

0.219 

0.261 

0  10  0 

7.461 

0.256 

7.200 

0.216 

0.264 

0  11  0 

7.717 

0.256 

7.420 

0.220 

0.266 

0  12  0 

7.988 

0.266 

7.641 

0.221 

Calculation  without  stretching  -  bending  coupling. 

Calculation  including  the  stretching  -  bending  coupling  (see  text). 
Empirical  model^^ 

Ab  initio  calculation^^ 


Figure  captions 


Fig.l  The  HOH  geometry. 

Fig.2  The  (ArEjO)"'’  adiabatic  PESs  dependence  on  the  HOH  valence  angle  for  different 
argon  positions  along  the  z-coordinate  (Fig.l). 

(a)  z  =  oo  (solid  lines)  and  z=4A  (dashed  lines). 

(b)  z=6A  (solid  lines)  and  z=5A  (dashed  lines). 

Fig.3  The  (ArH20)'*'  geometry. 

Fig.4  (a)  The  (ArH20)+  adiabatic  states  as  a  function  of  the  argon  position  p  (see 
Fig.3)  for  a  fixed  HOH  geometry  of  the  isolated  H2O  and  frozen  HOH  orientation, 
(b)  The  Ar  charge  along  p  (for  the  same  conditions  described  in  (a)). 

The  numbers  k=2,4,5  denote  the  Ar‘*‘H20  electronic  states  (see  text) 

Fig.5  The  Ar++H20  and  Ar+(H20)+  diabatic  energy  levels  for  fixed  HOH  geometry 
and  orientation.  Dashed  line  represents  the  Ar++H20  system  with  the  H2O  in 
the  electronic  —  vibrational  ground  state  (for  the  situation  given  in  Fig.3a).  Solid 
lines  are  for  the  various  excited  vibrational  bending  (OW2O)  states  of  the  A. 
(H20)+  ion  belonging  to  the  Ar+(H20)+  system  (for  the  situation  described  in 
Fig. 3b) .  The  results  are  for  impact  parameter  b=3A  and  for  the  following  geome¬ 
tries: 

(a)  The  H2O  atoms  coordinate  are:  0(0.065,0,0),  H(-0.520,0,±0.759)  and  the 
(H20)+  atoms  coordinates  are:  0(0,0,0),  H(0,0,±1.008). 

(b)  The  H2O  atoms  coordinate  are:  O (-0.065,0,0),  H(0.520,±0.759,0)  and  the 
(H20)+  atoms  coordinates  are:  0(0,0,0),  H(0,±1.008,0). 

(c)  The  H2O  atoms  coordinate  are:  0(0,0,-0.065),  H(±0.759,0,0.520)  and  the 

(H20)+  atoms  coordinates  are:  0(0,0,0),  H(±1.008,0,0). 


Fig.6  The  Ar+H20  geometry  for  the  full  scale  calculation 


Fig.7  (a)  Vibrational  bending  energy  levels  as  a  function  of  p  for  the  situation  given 
in  Fig. 6.  Dashed  line  is  for  the  electronic  state  correlated  asymptotically  with 
and  solid  lines  are  for  the  electronic  state  correlated  asymptotically 
with  Ar+(H20)+.  The  results  are  for  impact  parameter  b=3A.  (More  details  are 
given  in  the  text). 

(b)  The  valence  HOH  angle  and  Ar  atom  charge  as  a  function  of  p  for  the  sit¬ 
uation  given  in  Fig.6. 

—  •  —  •  HOH  valence  angle  (left  axis) 

-  Ar  charge  (right  axis) 
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